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1. Introduction 

This paper describes a computational method 
for the loop integrals. This method is named as 
Direct Computation Method and its advantage 
is handling singularities in a purely numerical 
way which appear in the denominator of the in- 
tegrand . We have already computed several 
diagrams with one-loop and two-loop such as 

• One-loop vertex and box diagrams with/ 
without infrared divergence [jl], ||, ^ Q], 

• Two-loop planar and non-planar vertex 
diagram [||]. 

In this paper, we present the brief description 
of this method and show another examples of 
loop integrals. 

2. Direct Computation Method 

This method is based on a numerical multi- 
dimensional integration to get the sequence of 
the integration approximations of the loop inte- 
gral, {/(£/)},/ = 0, 1,2, and an extrapolation 
technique for the convergence of the sequence. 
Here, the sequence of {£/},/ = 0, 1,2, ... is given 
as £/ = £o x a ~ l an d £o and a are real constants 
such as 256 and 2 respectively [jl]]. When we 
take the limit of £ — ► 0, we get the result of the 
loop integral. 

In Direct Computation Method, the inte- 
gration routine plays a central role. We have 
been using two different integration routines, 
DQAGE and DE. The former, DQAGE, is an adap- 
tive algorithm routine in QUAD PACK [jfj. We 
call the combination of DQAGE and the extrap- 
olation technique DQ-Direct Computation Method 
in this paper. The latter, Double Exponential 
formulas^, shortly DE, uses tanh(n/2*sinh(t)) 
transformation for numerical integration. We 
call the combination of DE and the extrapola- 
tion technique DE-Direct Computation Method 
in this paper. 



For both integration routines, we are com- 
bining Wynn's £ algorithm [JTj, |8]] or Aitken ex- 
trapolation to accelerate the convergence. 

3. Examples of the Computation 

3.1 One-loop box contributing to gg —> bbH 

The first example is the the one-loop box 
diagram with complex masses contributing to 
gg — ► bbH in Fig. [I] In the computation we 

Figure 1: Box diagram contributing to gg — > bi>H 

Graph 42 




produced by GRACEFIG 



takey^ = 353 GeV as an example and we fix 

= Jl^mf+M^) ps 271.06 GeV. As for 
the mass parameters, m t = 174 GeV, A% = 80.3766 
GeV and Me = 165 GeV. We introduce the com- 
plex masses as mj —> mj — im t T t and Mjy — ► 
M% - iM w F w with r, = 1.5 GeV and T w = 2. 1 
GeV. The results of the numerical computation 
agree perfectly to the analytic results reported 
by L. D. Ninh et al. [|] as in Fig. @ and Fig. | 

3.2 e + e~ — > ZZ — > e + e~Z 

The second example is the one-loop pen- 
tagon diagram in Fig. ^ The loop integral is 

r 1 5 1 

1= dx\dx2dxT ! dx^dx= ) 8{\ — V xA — 
Jo " £i £> j 

(3.1) 

Here, D(x) is given as 
D(x) =Y j m] 

£=1 
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Figure 2: Integration result of the real part with 
complex masses as a function of Js~2, Our results 
agree to the 5- or 6-digit accuracy of Ninh's. 



Figure 4: Diagram of e + e — ► ZZ — > e + e Z 
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Figure 3: Integration result of the imaginary part 
with complex masses as a function of y/s2. Our 
results agree to the 4-, 5- or 6-digit accuracy of 
Ninh's. 
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- X\X2S[ —X2X3S2 — — X4X5S4 

- X 5 XiS 5 -X l X 3 Si2-X 2 X4S 2 3 -X3JC5J34 

- X 4 XiS A5 - X 5 X 2 S 5 l , (3.2) 

where m\ = m 3 = m 4 = 0.5 x 10~ 3 GeV and 
m2=m= > = 91.00 GeV. As an example, one nu- 
merical result is shown in Table [l] with a set of 
parameters corresponding to one phase space 
point (Table ^). These agree to the results by 
T. Ueda et al. [[K]]. The elapsed time for the 
computation is about 9.3 hours for the real part 
using AMD Opteron 2.2 GHz CPU. 

3.3 Two-Loop Self-energy 

The third example is the two-loop self-energy 



Graph 2725 




produced by GRACEt-IG 



Table 1: Numerical result of one-loop pentagon. 





Result 


Error 


real 


0.4118519 x 10~ 13 


0.410 x 10" 17 


imag. 


-0.2336871 x 10~ 12 


0.619 x 10" 16 



Table 2: A set of parameters corresponding to one 
phase space point 



Parameter 


value 


5l 2 


100000.00000 


5l5 


-14146.0960752976 


S23 


-30471.3126018059 




32384.1496580698 


S45 


37833.5682283554 



diagram in Fig. |5| 



The formulae of the loop integral is [11] 



a 5 1 

(I6n 2 ) 2 l= / dx\dx2dxj,dx 4 dxs8{\ — Y[ %i) . 

Jo fr{ CD 

(3.3) 

Here, D(x) is given as 

D(x) = -p 2 (x 5 (xi +x 3 )(x 2 +x 4 ) 

+ (X\ +X 2 )x 3 X4 + (X3 +X 4 )xiX 2 ) 

+ CM 2 , (3.4) 
C = (Xi +X2+X3 +x 4 )x 5 + (xi +x 2 )(x 3 +X4), 
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Figure 5: Diagram of the two-loop self-energy 



-<Eh 



M 2 



i=i 



(3.5) 



We show the numerical results with two 
sets of the mass assignment shown in Table |3[ 
The results with the first set are compared with 



ones by Kreimer [12] and by Kurihara et al. 
fll3[ ] in Fig. |6| The results with the second set 



are compared with ones by Kurihara et al. [13], 



by Bauberger et al. [ 14] and by Passarino et al. 
JIB] ] in Fig. ^. In this computation, the elapsed 
time for the computation becomes longer around 
the singularities. 

Table 3: Set of mass assignment in GeV 



Set# 


in i 


rm 


ra 3 


m\ 


m 5 


1 


150.0 


150.0 


150.0 


150.0 


91.17 


2 


Vi.o 


V2.0 


V4.0 


75.0 


V3.0 



Figure 6: Numerical results of the loop inte- 
gral with Set#l mass assignment as a function of 
p 2 /150.0 2 . The results of [|l2| Q are divided by 
the value of p 2 . 



0.00035 
0.0003 
0.00025 
0.0002 
0.00015 
0.0001 
5e-05 


-5e-05 
-0.0001 



■ K Yuasa et ai. [real 




Y. Kurihara et al. [real 




D. Kreimer [real 




F. Yuasa et al. [imag 




Y. Kurihara [imag 




D. Kreimer [imag 






f± \ 















Figure 7: Numerical results of the loop integral 
with Set#2 mass assignment as a function of p 2 . 
The sign of our numerical results is conformed to 
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[re^l 






0.4 


Y. Kurihara et al. 


real 








S. Bauberger et al. 


real 






0.2 



, G. Passarino et al. 


real 






-0.2 


k 

- a. 








-0.4 
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imag 




-0.6 






Y.Kurihara et al. 


imag 










S. Bauberger et al. 


imag 




-0.8 
-1 




d 

i 


G. Passarino et al. 


imag 





4. Parallel Computation 

As we have described, the elapsed time for 
the computation becomes longer to get the se- 
quence {/(£/)} in a good accuracy when the 
singularity becomes steeper. To reduce the com- 
putation time, parallel computing technique is 
often used. We have developed the parallel code 
of Direct Computation Method. We used MP 1 1 
library which is widely used in the parallel com- 
puting environment. We evaluated the speedup 
of the parallel code of DQ - and DE- Compu- 
tation Method taking the integrals of two-loop 
vertex and one-loop box diagram as examples. 
The results of the speedup are shown in Ta- 
ble ||, Fig. H and Fig. ^. Both parallel codes 
show a good speedup behavior. 

5. Summary 

In this paper, we presented several numer- 
ical results of the loop integrals by Direct Com- 
putation Method. This method is based on the 
combination of the numerical integration and 
an extrapolation technique. To reduce the com- 
putation time we developed the parallel code 



Message Passing Interface 
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Table 4: Speedup of the parallel code of DQ-Direct 
Computation Method for the two-loop vertex dia- 
gram with singularities. We use PPC 970 and PC 
farm for the measurement. PPC 970 is IBM Blade 
CenterJS20 with PowerPC 970 2.2 GHz CPUs. PC 
farm consists of 16 Intel Xeon 3.06 GHz CPUs. 



# of CPUs 


PPC 970 


PC farm 


1 


1.00 


1.00 


2 


1.53 


1.50 


4 


2.99 


2.54 


8 


5.50 


4.54 


16 




13.31 



Figure 8: Speedup of the parallel code of DE- 
Direct Computation Method for the one-loop box 
integral with (DE-xdl-s-positive) and with- 
out (DE-xdl-s-negative) singularities. We 
use CRAY XD1 for the measurement. 




Figure 9: Speedup of the parallel code of DE- 
Direct Computation Method for the one-loop box 
integral with (DE-bg-s-positive) and with- 
out (DE-bg-s-negative) singularities. We use 
IBM BlueGene/L at KEK for the measurement. 




and the results of the speedup measurement are 
also shown. 
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